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ABSTRACT 

We have used the self-consistent vertical disc models of the solar neighbourhood pre- 
sented in Just & Jahreiss (2010), which are based on different star formation histories 
(SFR) and fit the local kinematics of main sequence stars equally well, to predict 
star counts towards the North Galactic Pole (NGP). We combined these four different 
models with the local main sequence in the filter system of the SDSS and predicted the 
star counts in the NGP field with b > 80 deg. All models fit the Hess diagrams in the 
F-K dwarf regime better than ±20 percent and the star number densities in the solar 
neighbourhood are consistent with the observed values. The x 2 analysis shows that 
model A is clearly preferred with systematic deviations of a few percent only. The SFR 
of model A is characterised by a maximum at an age of 10 Gyr and a decline by a factor 
of four to the present day value of 1.4M0/pc 2 /Gyr. The thick disc can be modelled 
very well by an old isothermal simple stellar population. The density profile can be 
approximated by a sech at function. We found a power law index a t = 1.16 and a scale 
height h t — 800 pc corresponding to a vertical velocity dispersion of cr t = 45.3 km/s. 
About 6 percent of the stars in the solar neighbourhood are thick disc stars. 

Key words: Galaxy: solar neighbourhood - Galaxy: disc - Galaxy: structure - 
Galaxy: evolution - Galaxy: stellar content - Galaxy: kinematics and dynamics 
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1 INTRODUCTION 

In Just & Jahreiss (2010) (hereafter Paper I) a self- 
consistent model of the vertical structure of the Milky Way 
disc in the solar neighbourhood was presented. The model 
is based on the star formation history (SFR), the age ve- 
locity dispersion relation (AVR), and a simple chemical en- 
richment model. The vertical density profiles of the stellar 
sub-populations are self-consistently calculated in the total 
gravitational potential of stars including the contribution 
of the gas component and the dark matter halo. The in- 
put parameters are selected and optimised to reproduce the 
velocity distribution functions fi(W) for the vertical veloc- 
ity component W along the main sequence (MS). It turned 
out that for each SFR the range of fitting AVRs is very 
small. On the other hand the SFR is not well determined by 
the local kinematics only. Therefore we discussed in Paper I 
four different SFRs with similar % 2 values to cover the range 
of possible functional shapes (models A, B, C, and D). The 
overall range of best fitting AVRs is still well restricted. Due 
to the different age distributions and consequently different 
relative contributions of the sub-populations as function of 
scale height the vertical density profiles of MS stars differ 
significantly at large distances z above the Galactic mid- 
plane. 



A combination of the local model with number densities 
of MS stars in the solar neighbourhood allows predictions of 
star counts at high Galactic latitude. Therefore we can use 
star counts of large surveys for an independent test of the 
local model. As a new additional result we expect to find 
restrictions on the SFR. As local normalisation along the 
MS we use Hipparcos stars complemented by the Catalogue 
of Nearby Stars (CNS4) at the faint end (JahreiB & Wie- 
len 1997). The model predictions will be compared to the 
North Galactic Pole (NGP) data of the Sloan Digital Sky 
Survey (SDSS), which has collected at present the largest 
and most homogeneous database comprising about 10 8 stel- 
lar objects in the Milky Way (Gunn et al. 1998; Abazajian 
et al. 2009). The SDSS photometry is in the ugriz filter 
system. Despite the enormous wealth of information that 
the SDSS database provides in the ugriz filter system, 
the majority of current observational and theoretical knowl- 
edge of resolved stellar populations is based largely on the 
Johnson-Kron-Cousins U BV Rclc photometric system and 
some other systems such as the Stromgren, DDO, Vilnius, 
and Geneva systems. To overcome this difficulty, we use an 
empirical transformation of the nearby MS star photometry 
into the ugriz system (Just & Jahreiss 2008). For each the- 
oretical model we derive a best fit solution for the full Hess 
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diagrams in (g — r, g), which quantify the number density 
distributions in the colour magnitude diagram (CMD). 

The paper is organised as follows. In Sect. 2 we de- 
scribe the data selection, in Sect. 3 important properties of 
the theoretical models are presented, in Sect. 4 the local 
normalisation is discussed, in Sect. 5 we present the best fit 
procedure, in Sect. 6 we discuss the results, and in Sect. 7 
we draw some conclusions. 



2 SDSS DATA 

The Sloan Digital Sky Survey (SDSS) is an imaging and 
spectroscopic survey (York et al. 2000) that has since 1998 
with its 2.5 m dedicated telescope mapped more than a 
quarter of the sky (Gunn et al. 1998). Photometric sky 
coverage of the SDSS Data Release 7 (DR7) amounts to 
ll,663deg 2 , including a 7, 646 deg 2 large contiguous area 
around the NGP (Adelman-McCarthy et al. 2008; Abaza- 
jian et al. 2009). 

Smith et al. (2002) defined the u g r' i z photometric 
system on 158 standard stars, a subset of UBVRcIc stan- 
dard stars from Landolt (1992), using the USNO-l.Om tele- 
scope. Unfortunately, the photometric system of the 2.5m 
SDSS telescope, denoted as ugriz, slightly differs from the 
u' g' V i' z' one (Abazajian et al. 2003). Moreover, the SDSS 
standards are too bright and saturate in the 2.5m SDSS tele- 
scope during its normal operational mode. These inconve- 
niences have been resolved by using fainter secondary stan- 
dards scattered throughout the SDSS survey area and by 
using simple linear transformation equations between the 
two photometric sets (Tucker et al. 2006). The nightly pho- 
tometry obtained by the SDSS 2.5m telescope can be thus 
calibrated to the native ugriz system with magnitude zero 
points accurate on the AB system to within a few percent. 

Imaging data are produced simultaneously in the five 
photometric bands, namely u, g, r, z, and z (Fukugita et al. 
1996; Gunn et al. 2006; Smith et al. 2002; Hogg et al. 2001). 
The images are automatically processed through specialised 
pipelines (Lupton, Gunn, & Szalay 1999; Lupton et al. 2001; 
Stoughton et al. 2002; Pier et al. 2003; Tucker et al. 2006) 
producing corrected images, object catalogues, astrometric 
solutions, calibrated fluxes and many other data products. 
SDSS photometry is homogeneous and deep (r < 22.5), re- 
peatable to 0.02 mag (ivezic et al. 2003) and with a zero- 
point uncertainty of ~ 0.01 — 0.02 (Abazajian et al. 2004; 
Ivezic et al. 2004). In DR7 a homogeneous photometry over 
the full sky coverage was established at the 1 percent level 
in griz and 2 percent in u in a process called ubercalibration 
(Padmanabhan, N., et al. 2008). 

We have limited our analysis to the stars in the mag- 
nitude range 14 ^ g ^ 20. 5 1 . For brighter stars the CCD 
camera of the SDSS telescope saturates. At the faint end 
we wanted to completely avoid the problems in the galaxy- 
star separation and we set a conservative magnitude limit 
for this purpose. The SDSS photometric data contain also 
quality flags for each object to aid in the selection of "good" 

1 Magnitude limits applied to the dereddenend magnitudes. We 
applied the standard dereddening of the data based on the ex- 
tinction map of Schlegel et al. (1998) as given in the data base of 
SDSS. 



measurements 2 . We have carefully analysed the appearance 
of these photometric flags with respect to the object bright- 
ness. We concluded that the problematic flags relate mostly 
to the stars fainter than our faint magnitude limit. Many 
of these flags are also tightly correlated with the magnitude 
measurement error. To make our stellar samples as com- 
plete as possible we have finally applied only one " cleaning" 
criterion, i.e. we have rejected all the measurements with 
reported magnitude error larger than 0.2 in g or in r filter, 
which typically accounted for 1 percent or even less of the 
total star counts. It is also necessary to mention that the 
typical magnitude error in the chosen magnitude range is 
much smaller than the applied magnitude error limit. 

We selected the colour range —0.2 < g — r < 1.2 with a 
symmetric bin size of ±0.025 mag. This range covers the MS 
down to K dwarfs, where the local normalisation is reliable. 
In order to test the predictions from the vertical structure 
of the local disc model (Paper I) with the available SDSS 
photometric data we have selected a field at the NGP. The 
field should not reach too low galactic latitude in order to 
avoid projection effects and a dependence of the star counts 
on the radial properties of the disc model. At the same time 
the field should be as large as possible in order to reduce 
the Poisson noise of the star counts. We have chosen a NGP 
field with Galactic latitudes b > 80° and an area of Angp = 
313.36 deg 2 . In g magnitude we use a resolution of dg = 
0.01 mag, but for the fitting procedure we applied a car-box 
smoothing of Ag =0.1, 0.2, and 0.5 mag, where the last value 
is our standard case. The resolution in colour is d(g — r) = 
0.05 mag. The Hess diagram of the data is shown in the 
bottom panel of figure 4. 

In order to avoid confusion of observational errors in 
the SEGUE data and uncertainties of the model predictions 
we compare dereddened number counts Y{g — r, g) in the 
colour - apparent magnitude plane. 



3 THE MILKY WAY MODEL 

We use the self-consistent local model of the thin and thick 
disc described in detail in Paper I and add a simple stel- 
lar halo component. In this section we report the necessary 
features to understand the procedure applied here. The disc 
model relies on the kinematics of MS stars in the solar neigh- 
bourhood measured by volume complete local samples of 
stars combining Hipparcos stars and faint stars of the CNS4. 
The properties of the model depend on the total stellar mass 
density in the solar neighbourhood and only weakly on the 
adopted Initial Mass Function (IMF). It is independent of 
star counts at large distances. 

The thin disc model at the solar radius is based on a 
pair (SFR,AVR) as input functions, which determines the 
age distribution and the velocity distribution functions of 
the stellar sub-populations. The vertical density profiles are 
calculated self-consistently in dynamical equilibrium in the 
total gravitational potential of stars, gas and dark matter 
halo. The chemical enrichment is also included for the de- 
termination of MS lifetimes, colour indices and luminosities 

2 See the following link for the clean photometry prescriptions: 
http: / / cas.sdss.org/ dr7/en/help/docs/realquery.asp#flags. 
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model A: g-r > 0.40 (>12. Gyr) 
= 0.35 (= 6.0Gyr) 
= 0.30 (= 3.0Gyr) 
= 0.25 (= 1.6Gyr) 



age [Gyr] 

Figure 1. The upper panel shows the SFRs as function of time 
and the lower panel gives the AVRs of models A - D as function 
of age (age is running backwards). 



from Padova population synthesis models. For each pair of 
functions (SFR,AVR) the local disc model provides a unique 
connection of local star counts and the number density of 
stars above the galactic plane with no additional free pa- 
rameters. In Paper I we compared four different models A, 
B, C, D which all yield similar best fit \ 2 values. The SFR 
and AVR of these models are plotted in figure 1 . The density 
profile of MS stars at a given colour index is characterised 
by the MS lifetime of the stars. It is composed of a series 
of density profiles of coeval stars according to the SFR up 
to the lifetime. The top panel of figure 2 shows the nor- 
malised density profiles of model A for the relevant colour 
index range in g — r. Bluer thin disc stars are too bright to 
be visible in SDSS at high Galactic latitude. 

The lower panels of figure 2 show the density profiles of 
models A - D divided by a sech 2 (z/2z s ) profile with expo- 
nential scale height z s for two colour (MS lifetime) bins. For 
both lifetimes we use the corresponding exponential scale 
height of model A, i.e. z s =270, 108 pc respectively. The mid- 
dle panel is in linear scale, whereas the lower panel for stars 
with shorter lifetime is in log-scale. The deviations from a 
sech 2 -profile of model A are significant for all lifetimes, be- 
cause the populations are not isothermal and their spatial 
distribution is influenced by the gravitational potential of all 
other components. For younger populations the deviations 
are much stronger than for older populations. Between the 



q 0.6 




z [kpc] 

Figure 2. The top panel shows normalised density profiles of 
model A for different colour bins (with the corresponding MS 
lifetime in parenthesis). The middle panel shows the deviations 
of the density profiles for models A - D from a posech 2 (z / 2 z s ) 
profile with z s = 270 pc for stars with lifetime larger than 12 Gyr 
(g — r > 0.4). The lower panel is in log-scale and shows the same 
for a lifetime of 3 Gyr using z s = 108 pc. 



different models the differences exceed a factor of two. Mod- 
els B-D with a lower fraction of stars older than 10 Gyr 
require larger velocity dispersions of the old thin disc popu- 
lations (see figure 1). This leads to shallower density profiles 
at z > 1 kpc. All profiles are based on a Scalo IMF. The 
differences at the mid-plane z = can be corrected by ad- 
justing the IMF, which will be done implicitly by fitting the 
local normalisation (see below) . In any case large differences 



4 A. Just et al. 



of star counts as function of distance (apparent magnitude) 
remain. 

In each g — r colour bin the density profile is charac- 
terised by the MS lifetime. According to Paper I and Just 
& Jahreiss (2008) we use for the bins g - r=0.15, 0.2, 0.25 
the lifetime 1.2, 1.4, 1.6 Gyr, respectively. In the colour bins 
g — r=0.3, 0.35, 0.4 we include the contribution of turnoff 
stars and use lifetime ranges 2.4-3.2, 4.0-6.2 and 8.2-12 Gyr, 
respectively. For all colours g — r ^0.45 the lifetime is larger 
than 12 Gyr leading to the same density profile. 

An isothermal thick disc component is self-consistently 
included in the local model. We have shown in Paper I that 
it can be well fitted by a 



pt(z) = p t ,osech at [z/(a t z t )] 



(i) 



law with local density p t ,o and exponential scale height z t . 
Since the influence of the thick disc on the total gravitational 
potential &(z) is very small, the thick disc parameters can 
be varied in a large range with negligible influence on &(z) 
and the thin disc structure. From the Jeans equation it fol- 
lows that zt and a^ 1 are both proportional to the square of 
the velocity dispersion a t of the thick disc in a given gravi- 
tational potential. As a consequence we find 



atZt ~ const. 



(2) 



which we use to maintain dynamical equilibrium, when 
changing the thick disc parameters. For the stellar popu- 
lation of the thick disc we adopt a simple population with 
an age of 12 Gyr and metallicity [Fe/H]=-0.7. 

We include a simple stellar halo described by a flattened 
power law distribution 



Ph — Ph,o 



R 2 + z 2 /b 2 
R 2 



(3) 



with flattening b and local normalisation ph,o- For the dis- 
tance of the Sun to the Galactic centre we adopt Ro = 
8.0 kpc. Since we are looking only to one line of sight and to 
distances not large compared to Ro, the parameters and 
b are strongly degenerated. A stronger flattening requires a 
shallower slope in order to reproduce similar stellar densities 
at a distance of z = 5 — 10 kpc. One choice to get minimum 
X 2 values is ah = —3.0 and b — 0.7, which we fix for the 
further investigations. 

For the calculation of star counts we use the density 
profiles of each component v at distance s along the line 
of sight pointing to the Galactic coordinate position (l,b) 
normalised at the solar position (R,z)q = (Rq,0) 



Po 



with p = p(Ro, 0), 



(4) 



z — s sin ( 



R — y / ' Rq — 2Ros cos b cos I + s 2 cos 2 b. (5) 

For the thin disc the density profile depends also on colour 
g — r. The vertical offset zo ~ 20 pc of the solar position 
can be easily included but is negligible for high Galactic lat- 
itude fields. The density distributions of thin and thick disc 
can be extended in the radial direction by an exponential 
profile, if necessary. The normalised density profiles of each 
component are transformed to number density profiles by 
multiplying with the local number densities n u ,o(g — r, M g ) 
determined from the HRD in the solar neighbourhood. 
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Figure 3. The local MS of thin disc (red full circles and line), 
thick disc (green open circles), and stellar halo (blue triangles). 
The MS stars of the CNS4 and the fiducial sequence of M13 with 
metallicity [Fe/H]=-1.53 are over-plotted. 



4 LOCAL NORMALISATION 

The Galaxy model described in the last section must be 
complemented by the local stellar number density n u ,o(g — 
r, M g ) of each component v. In general the full HRD of a 
volume complete sample should be used. Since most of the 
modelled CMD is dominated to more than 95 percent by 
MS stars, we restrict the star count predictions to MS stars 
including a correction for turnoff stars. Only at the bright 
red corner of the CMD the K giants of thick disc and halo 
contribute significantly. 

The stellar content in the solar neighbourhood is quan- 
tified by the number of MS stars N25 with 



N 25 = J2 N25 >»( g ' 



with v — s,t,h (6) 



in a sphere of 25 pc radius with volume V25 = 65,450pc 3 . 
Here the contributions of thin disc, thick disc, and halo are 
denoted by the indices s, t, and h, respectively. In Just & 
Jahreiss (2008) absolute magnitudes and colours of the mean 
MS in the u, g, r, i, z filter system were determined based on 
different transformation formulae available in the literature. 
We use that mean MS for the thin disc, because the influence 
of thick disc and halo stars in the solar neighbourhood is 
negligible. 

In figure 3 we show M g (g — r) for the thin disc based on 
the transformation J06 (Jordi et al. 2006). The black crosses 
denote the MS stars of the CNS4 which are the basis for the 
determination of the mean MS. Each red full circle along 
the MS denotes a data point, where N25, s (g — r,M E ) is a 
free fitting parameter. The locus of the MS corresponds to 
a slightly sub-solar mean metallicity and is consistent with 
the fiducial sequences of the open clusters M67 ([Fe/H]=0) 
and NGC2420 ([Fe/H]=-0.37) of An et al. (2008). Since the 
mean metallicity of the thin disc population decreases with 
increasing distance from the mid-plane, a correction to the 
mean luminosity in each colour bin is necessary. We use a 
simple analytic approximation 



AM g ([Fe/H](z)) = 0A 



z + 600pc 



mag, 



(7) 
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which is consistent with the determination of Ivezic et al. 
(2008) for the M r luminosity of thin disc MS stars. 

For the thick disc we use an interpolated MS with 
N 2 5,t(g — r,M g ) corresponding to an old population with 
metallicity [Fe/H] « —0.7 (green open circles in figure 3). 

The stellar halo is represented by N 2 5^(g — r, M g ) corre- 
sponding to an old population with metallicity [Fe/H]= —1.5 
(blue triangles in figure 3), where we used the fiducial se- 
quence of M15 (An et al. 2008) with an extrapolation to the 
faint end. 

In the turnoff regime of F stars we split the MS lumi- 
nosity into two or three points with different luminosities in 
order to include the brighter turnoff stars. The total number 
of thin disc, thick disc and halo should add up to the ob- 
served total number of MS stars N 2 $(g — r), which were also 
determined in Just & Jahreiss (2008) (histogram in figure 
8). We treat the local number densities of the components 
as free fitting parameters and then compare the results of 
the best-fitting model to the observed stellar content in the 
solar neighbourhood. 

The predicted Hess diagrams, i.e. the star counts in the 
CMD N(g — r,g) per d(g — r)dg bin, in a cone with cross 
section dl cos bdb pointing to the Galactic coordinate posi- 
tion (l,b) are calculated by adding up the contributions of 
each component along the line of sight. The contribution of 
each component v by a volume element at distance s, using 
ds/s = 0.21og(10) dg, is given by 

AN u dg = N 2 5,u(g - r,M g )p„(s)s 2 dsdl cos bdb (8) 
N 2 5,v(9 ~ r, M g ) Piy (s)s 3 dg 4tt 

V 25 2.17 41,253deg 2 

g = M g + AM g + 5 log 10 s - 5. 

Adding up the contributions of all components along the line 
of sight yields the predicted number N(g — r, g) in each bin of 
the Hess diagram. In the next step the values are smoothed 
in the same way as the SDSS data. The star numbers in the 
Hess diagrams are normalised to 1 deg 2 at the sky (second 
line of equation 8). In (g — r,g) we need to distinguish be- 
tween the resolution d(g — r)dg, the smoothing Ag and the 
normalisation. All Hess diagrams are normalised to 0.1 mag 
in (g — r) and 1.0 mag in g. 



5 FITTING PROCEDURE 

The best-fit procedure in the colour bins are independent 
of each other. In each colour bin the contribution of thin 
disc, thick disc, and stellar halo to the local star counts in 
V25 are quantified by the free parameters N 2 5,s, A^5,t and 
A^25,h, respectively. The (g — r) arguments are dropped here 
for simplicity. We use the nonlinear Levenberg-Marquard 
algorithm (Press et al. 1992) to minimise in each colour bin 
i 

2 (logjjj -logA^ij) 2 ... 2 1 , Q x 

Xi = 2^ —2 wlth ^ ( 9 ) 

in log-scale with Poisson noise for the statistical weights 
(Jij . In linear scale the regions with maximum density would 
dominate the % 2 value resulting in an unsatisfactory overall 
fit. The x? are derived from normalised values and have to be 
corrected by the area Angp in units of deg 2 and the bin size 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 4. Top to bottom: Contributions of thin disc, thick disc 
stellar halo, and full Hess diagram of model A-R06-05 (see table 1) 
followed by the NGP data. Colour coding ranges from 1 (purple) 
to 100/deg 2 (red) in log-scale. 



d(g — r)Ag in units of 0.1 magx 1 mag. Due to the smoothing 
in g over n data points xf l n is the average over all subsets 
of independent bins with % 2 k shifted by n x dp. In each 
colour bin the degrees of freedom dof\ for the independent 
bin sets depend on the fit regime in g, the bin size Ag and the 
number of fitting parameters. We derive the mean reduced 
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Table 1. Parameters of the best fit solutions. 



2 _ Angp d(g 



■r) Ag 



1 deg 2 0.1 mag 1 mag n x dof\ 



Xi 



(10) 



In order to get reliable fits we restrict the colour range of 
the fitting regime. For the thin disc we fix N25,s(g — r, M g ) 
for g — r < 0.35, since the main contribution falls outside 
the bright limit g — 14 mag. For the halo we extrapolate 
N25,h(g — r, M g ) for g — r > 1.0, since the main contribution 
falls outside the faint limit g — 20.5 mag. Additionally the 
part of the CMD, which we use to minimise \ 2 is restricted 
dependent on the aspect of investigation. 



6 RESULTS 

We discuss first the construction of the Hess diagrams and 
compare the best fit results of models A - D. Then we in- 
vestigate the dependence of the fitting result on the fitting 
regime and the MS properties for model A. 



6.1 Hess diagram fitting 

The bottom panel of figure 4 shows the Hess diagram 
of the NGP data with a total of 276,180 stellar objects. 
Star number densities normalised to 0.1 magx 1.0 mag in 
(g — r, g) and 1 deg 2 are colour coded in log-scale ranging 
from 1. . . 100 (purple ...red). The plots are smoothed in 
colour and magnitude in boxes in steps of 0.01 mag with box 
size A(g — r) x Ag =0.05 magx 0.5 mag (in the fitting pro- 
cedure the colour bins are not smoothed). Most of the Hess 
diagram is strongly dominated by MS stars. Exceptions are 
turn-off stars in the colour range g — r w 0.3-0.4 and red gi- 
ants of thick disc and/or halo in the upper right corner. The 
nature of the faint very blue objects with g — r <0.15 is un- 
clear (miss-identified extragalactic sources, White Dwarfs, 
halo BHB stars, . . . ). 

The three top panels of figure 4 show the contributions 
of the different components to the full Hess diagram for 
model A. There is a significant overlap of thin/thick disc, 
as well as thick disc/halo. Therefore the stellar halo must 
be included to determine the thick disc properties and the 
thick disc is needed to fix the thin disc parameters. The sec- 
ond last plot shows the sum of all components for model 
A. 

For the comparison of disc models A - D we restrict 
the colour range to 0.5 ^ g — r ^ 1.2 mag in order to avoid 
the F turn-off regime, which may be improperly modelled 
and thus dominate the total \ 2 • The magnitude range is 
14.5 ^ g ^ 20.5 mag for safely excluding a contribution of 
saturated stars brighter than 14.25 mag in the brightest bins 
and a significant contamination by misidentified extragalac- 
tic sources at the faint end. In all models the thick disc 
parameters are varied according to equation 2 to find the 
minimum % 2 value. The parameters of the best fits are col- 
lected in the first four rows of table 1. In the model names 
there are tokens for the local normalisation (-J06 for Jordi 
et al. (2006), -R06 for Rodgers et al. (2006), -C08 for Ono- 
nis & Gaskell (2008)) and the smoothing Ag (-01, -02, -05, 
respectively) attached. For the first block '-r' is attached for 
the reduced fit regime. 



model 


Ag[mag] 


x 2 


^25 


zt [pc] 


cut 


at [km/s] 


A-J06-05-r 


0.5 


2.59 


739 


800 


1.16 


45.3 


B-J06-05-r 


0.5 


3.75 


748 


880 


1.07 


47.4 


C-J06-05-r 


0.5 


3.84 


811 


885 


1.23 


51.2 


D-J06-05-r 


0.5 


5.64 


709 


930 


0.99 


48.0 


A-R06-05-r 


0.5 


2.61 


715 


800 


1.16 


45.3 


A-C08-05-r 


0.5 


2.60 


733 


800 


1.16 


45.3 


A-J06-05 


0.5 


4.50 


762 


800 


1.16 


45.3 


A-R06-05 


0.5 


4.31 


710 


800 


1.16 


45.3 


A-C08-05 


0.5 


5.09 


752 


800 


1.16 


45.3 


A-C08-02 


0.2 


2.53 


750 


800 


1.16 


45.3 


A-C08-01 


0.1 


1.84 


753 


800 


1.16 


45.3 



Note. The models are named by the disc model A - D followed by 
the local normalisation and the smoothing Ag. For the first block 
V is added for the reduced fit regime. The observed number of 
stars range from N 25 = 726 (with R06) to N 25 = 770 (with C08). 



6.2 The star formation history 

The thin disc models A - D have very different SFRs leading 
to different star count predictions. In each model we opti- 
mised the thick disc scale height and power law index ac- 
cording to equation 2 starting from a self-consistent isother- 
mal model. In figure 5 the Hess diagrams of the models 
are presented in the left column (A - D: top - bottom). 
The middle column shows the relative differences (data- 
model) /model, where the coloured region covers the range 
-20 percent (purple) to +20 percent (red) in linear scale. 
Lower values are represented by black colour and higher 
values by white colour. The right column shows the con- 
tribution of each bin to x? (equation 9) in log-scale. 

All models fit most of the Hess diagram better than 
±20 percent. The local normalisation of the blue part with 
g — r < 0.5, which is not included in the best fitting, is 
adapted by hand to model A and not varied for the other 
models. In the bright red triangular region there is a sig- 
nificant fraction of stars missing in the models as expected, 
since the SDSS data contain red giants of thick disc and 
halo, which cover that colour-magnitude regime. 

In figure 6 the contributions of thin disc, thick disc, and 
stellar halo for model A are shown in the colour bin g — r = 
0.7 as a typical example. The maximum contribution of the 
thin disc is at g = 16.5 mag corresponding to z — 900 pc and 
of the thick disc at g = 19 mag corresponding to z = 2500 pc. 
These large distances demonstrate that it is very important 
to construct models, which hold also for the outer profiles 
of the components. 

Model A-J06-05-r is the best model with a reduced 
X 2 = 2.59 (see table 1). In model A the deviations of data 
and model are less than ±10 percent over the full colour 
range and in the magnitude regime, where the MS stars of 
thin disc, thick disc or halo dominate. The x 2 distribution 
in figure 5 shows a strong noise component and the con- 
tribution by systematic deviations of the predicted density 
profiles from the real profiles. In the top panel of figure 7 
the reduced xf values in each colour bin are quantified. The 
contributions to the total x? are fairly uniform. 

The main feature of the systematic discrepancies is the 
shallow valley in the transition of thin disc and thick disc 
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Figure 5. Comparison of models A,B,C,D-J06-05-r (top to bottom). From left to right: Hess diagram (same colour coding as in figure 
4), relative difference to SDSS data ranging from -0.2 (purple) to +0.2 (red) in linear scale, xfj distribution in log-scale. 



in the data, which is not reproduced by the model. This is 
a sign of an even steeper slope of the thin disc above 1 kpc. 
A continuous transition between the old thin disc with a 
maximum vertical velocity dispersion of a e — 25km/s and 
the thick disc with a t = 45.3 km/s can be excluded, because 
it would add more stars to the transition regime. 

In models B, C, and D, the fits in the magnitude range 
of dominating thin disc contribution are significantly worse, 
whereas the thick disc parameters can be adjusted to reach 
a similar good fit at fainter magnitudes as in model A. This 
is also quantified in the % 2 values (see table 1). Model D 
with the small disc age and strongly declining SFR is the 
worst model with x 2 = 5.64. Models B and C are compara- 
ble with x 2 — 3.75 and 3.84, respectively. All three models 
overestimate the thin disc density at large heights leading 
to a deficit closer to the plane. A more detailed look on the 
Xi distributions shows that model A is superior in a wide 
colour range (figure 7). 

Since we used the local number densities of all compo- 
nents as free fitting parameters, a comparison of the fitting 



results and the observed star counts in the solar neighbour- 
hood is a crucial test of the models. In figure 8 the his- 
togram with error-bars shows the data from the CNS4. The 
full (blue) circles are the results of model A - D (from top 
to bottom) . In all cases there is a reasonable match and also 
the total number of MS stars agree (see table 1). The ob- 
served number of MS stars range from N25 = 726 (with R06) 
to N 2 5 = 770 (with C08). Additionally the contribution of 
thin disc, thick disc and halo is plotted separately in figure 
8 (enhanced by some factor to make it visible) demonstrat- 
ing that the local number densities of all components are 
smooth functions of (g — r). A more detailed discussion is 
given in section 6.3. 

As a bottom line we find that the local disc model A 
is fully consistent with the NGP star count data of SDSS. 
Model A with a relatively large fraction of stars older than 
8 Gyr and a correspondingly small maximum vertical veloc- 
ity dispersion of a e = 25 km/s is the best fitting model. In 
the next sub-sections we test the robustness of the result 
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shows the comparison of disc models A - D. The bottom panel 
shows the dependence on the filter transformation and on the 
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Figure 8. Local normalisations for disc models A - D (top to 
bottom). The histogram shows the data from the CNS4 and the 
full blue circles are the best fit values. The contributions of thin 
disc, thick disc (x5), and stellar halo (x50) are shown separately. 



with respect to the filter transformations used for the solar 
neighbourhood and with respect to the smoothing in g. 



6.3 Filter transformations 

In figure 9 the local MS is shown for three different transfor- 
mations (J06 for Jordi et al. (2006), R06 for Rodgers et al. 
(2006), C08 for Chonis & Gaskell (2008)). There are small 
systematic deviations in M g (g — r) which may influence the 
best fit parameters of the model. The histograms in figure 
11 show the local star counts for R06 and C08. Additionally 
to the differences in the absolute magnitudes for thin and 
thick disc the fraction of thin disc turn-off stars is reduced in 
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Figure 9. Local CNS4 stars and mean MS using different fil- 
ter transformations. Each star is plotted three times in different 
colours to show the systematic shift of the MS. 

C08 and is set to zero in R06 in order to test their influence 
on the x 2 values and number counts. 

We present the comparison using the different local MS 
for model A. The effects in models B, C and D are very 
similar. First we derived the best fit models for the reduced 
colour-magnitude regime as before in A-J06-05-r for models 
A-R06-05-r and A-C08-05-r (see table 1). The total % 2 values 
and local normalisations of the three models are very similar. 

Then we extended the fit regime to the full CMD and re- 
calculated the local normalisations and the x 2 values. Figure 
10 shows the relative deviations of data and models. The 
patterns of the deviations differ slightly in the regime, where 
the thin disc dominates. The corresponding x? values in the 
colour bins are shown in the lower panel of figure 7. For 
the extended fit regime there are small differences arising 
mainly from the turn-off stars. In the colour bin g — r — 0.4 
the fraction of turn-off stars is 80 percent (J06), percent 
(R06) and 60 percent (C08), respectively. Model A-R06-05 
without turn-off stars gives the best result. A consequence is 
that a larger fraction of stars are allocated to the thick disc 
instead of the thin disc in model A-R06-05 (see figure 11). At 
g — r = 0.35 model A-J06-05 with 55 percent turn-off stars 
gives the best result and for g — r ^ 0.3 the turn-off stars 
are too bright to contribute to the fit. At the red end xf is 
influenced by the missing giants. Therefore it is not useful 
to discuss here the differences of the fits in more detail. 

The local normalisations do not differ significantly (see 
figure 11). The differences in the histograms for the transfor- 
mations R06 and C08 show the possible variations in the lo- 
cal star counts due to noise and the effect of a different slope 
in the MS leading to a shift of stars in colour. All transfor- 
mations are consistent with the local star count data. Only 
in the transition of K to M dwarfs the large number of thin 
disc stars in g — r — 1.2 may be a hint that the MS turning 
point is relatively blue as in C08. 

6.4 Smoothing 

Smoothing of the data has two main effects. On one hand 
statistical noise is reduced. On the other hand physical fea- 
tures in the data are smeared out and may be shifted sys- 
tematically. Figure 12 shows the NGP data and correspond- 
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Figure 10. Relative differences of NGP data and models A- JOG- 
OS, A-R06-05, and A-C08-05 with different filter transformations 
(top to bottom). 

ing models A-C08-05,02,01 for different smoothing lengths 
Ag =0.5, 0.2, 0.1 mag, respectively. The x 2 values are de- 
creasing due to to the increasing degrees of freedom (see 
table 1 and figure 7). The local normalisations are statisti- 
cally indistinguishable. Only in the colour bin g — r — 0.4 
a few percent of thin disc stars are shifted to the thick disc 
at higher resolution reducing the contribution to xf by the 
bright magnitude bins dramatically. 



7 CONCLUSIONS 

We have used four different models with different SFRs of 
the thin disc, which fit the local kinematics of main sequence 
(MS) stars (Paper I) and compared the star count predic- 
tions for the North Galactic Pole (NGP) field with b > 80° 
of the Sloan Digital Sky Survey (SDSS). For the thin disc we 
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Figure 11. Local normalisations for disc model A using different 
filter transformations. 



applied the absolute magnitudes of the local MS as deter- 
mined in Just & Jahreiss (2008). A self- consistent isother- 
mal thick disc and a simple stellar halo model was added to 
complete the contributions of MS stars to the star counts. 
We used the local normalisations of thin disc, thick disc and 
halo in each colour bin as free parameters and minimise x 2 
in the Hess diagram over a large colour-magnitude range in 
(9-r,g). 

All models match the observed number densities better 
than ±20 percent proving the reliability of the thin disc den- 
sity profiles by the self- consistent disc models in the range of 
\z\ < lkpc. The derived local normalisations are consistent 
with the star count data of the CNS4 in the solar neigh- 
bourhood. The x 2 analysis shows that model A is clearly 
preferred with systematic deviations of a few percent only. 
The SFR of model A is characterised by a maximum at an 
age of 10 Gyr and a decline by a factor of four to the present 
day value of 1.4 MQ/pc 2 /Gyr. In the thin disc the present 
day fraction of stars older than 8 Gyr is with 54 percent sig- 
nificantly higher than for models B, C, and D (Paper I). 
Especially model C with a constant SFR and model D with 
a disc age of 10 Gyr can be ruled out. 

The thick disc can be modelled very well by an isother- 
mal simple stellar population. The density profile can be 
approximated by a sech at (z / atht) function. For model A 
we find a power law index of a t = 1.16 in-between an expo- 
nential profile and a sech 2 -profile (the latter corresponds to 
an isolated isothermal disc). The exponential scale height is 
h t — 800 pc corresponding to a vertical velocity dispersion 
of at = 45.3 km/s. About 6 percent of the stars in the so- 
lar neighbourhood are thick disc stars. In Juric et al. (2008) 
the stellar density distribution in the Milky Way based on 
SDSS star counts was fitted by exponential thin and thick 
disc profiles. The result is a much larger thick disc scale 
height, which balances the flattening of the profile at low z. 

The results do not depend significantly on the filter 
transformations used for the local MS nor on the smoothing 
of the data in luminosity. For the future an extension of the 
model to include turn-off stars in more detail and the con- 
tribution of giants as well as a higher resolution in colour at 
the blue end would be very useful. 

We also plan to apply the model to lower Galactic lati- 
tudes in order to determine radial scale lengths of thin and 



thick disc as well as radial gradients in the stellar popula- 
tions. 
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